arXiv:1503.08865vl [astro-ph.GA] 30 Mar 2015 


Mon. Not. R. Astron. Soc. 000, [l HTo| (2014) Printed 1 April 2015 (MN lAT^X style file v2.2) 


Megahertz peaked-spectrum sources in the Bootes field I - 
a route towards finding high-redshift AGN? 

Rocco Coppejans 1 *, David Cseh 1 , Wendy L. Williams 2,3 , Sjoert van Velzen 4 and 
Heino Falcke 1,3,5 

1 Department of Astrophysics/IMAPP, Radboud University Nijmegen, P.O. Box 9010, 6500 GL Nijmegen, The Netherlands 
2 Leiden Observatory, Leiden University, PO Box 9513, 2300 RA, Leiden, The Netherlands 
3 Netherlands Institute for Radio Astronomy (ASTRON), PO Box 2, 7990 A A, Dwingeloo, The Netherlands 

4 Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, MD 21218, USA 

5 Max-Planck-Institute fur Radioastronomie, auf dem Hiigel 69, 53121 Bonn, Germany 


ABSTRACT 

We present a 324.5 MHz image of the NOAO Bootes field that was made using Very 
Large Array (VLA) P-band observations. The image has a resolution of 5.6 x 5.1 arcsec, 
a radius of 2.05° and a central noise of ~ 0.2 mJy beam -1 . Both the resolution and 
noise of the image are an order of magnitude better than what was previously avail¬ 
able at this frequency and will serve as a valuable addition to the already extensive 
multi wavelength data that are available for this field. The final source catalogue con¬ 
tains 1370 sources and has a median 325 to 1400 MHz spectral index of -0.72. Using a 
radio colour-colour diagram of the unresolved sources in our catalogue, we identify 33 
megahertz peaked-spectrum (MPS) sources. Based on the turnover frequency linear 
size relation for the gigahertz peaked-spectrum (GPS) and compact steep-spectrum 
(CSS) sources, we expect that the MPS sources that are compact on scales of tens 
of milliarcseconds should be young radio loud active galactic nuclei at high (z > 2) 
redshifts. Of the 33 MPS sources, we were able to determine redshifts for 24, with an 
average redshift of 1.3. Given that five of the sources are at z > 2, that the four faint 
sources for which we could not find redshifts are likely at even higher redshifts and 
that we could only select sources that are compact on a scale of ~ 5 arcsec, there is 
encouraging evidence that the MPS method can be used to search for high-redshift 
sources. 
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1 INTRODUCTION 

The National Optical Astronomy Observatory (NOAO) 
deep wide-held survey Bootes held covers ~ 9 deg 2 in the 
optical and near infrared Bw, R, I and K bands (Jannuzi 


|fe Dey|1999[). In addition, the held has also been covered in 


2004 


the X-ray (Murray et al. 


2005 Renter et al. 20051, ultra¬ 


violet (Martin et al. 20031, mid infrared (|Eisenhardt et al. 


|Martin et al.||2003|) and radio: 3.1 GHz (|Croft et al. 


20131, 1.4 GHz (de Vries et al. 
153MHz (Williams et al.| |2013 


2002 Higdon et al. 20051, 


Interna et al. 2011) and at 


62 MHz, 46 MHz and 34 MHz by |van Weeren et al. | ( f2014| ) . 
Furthermore, the active galactic nuclei (AGN) and Galaxy 
Evolution Survey (AGES; Kochanek et al. "2012) have mea¬ 
sured redshifts for 23,745 galaxies and AGN in the held. 
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This extensive multiwavelength data set have in the past 
been used to study, amongst other things: the properties 

the spacial clustering 


of galaxy groups (Vajgel et al. 


of quasars (Hickox et al. 20111, search for radio transients 
(Bower etTa l.J2011R study the start formation rate in optical 


20141 


obscured galaxies ( Calanog et al.|2013 l and study the com¬ 
position of gas and dust in AGN ( Usman et al.|2012 |. The 
325 MHz observations of the held presented here will serve 
as a valuable addition to the multiwavelength coverage of 
the held and will, as shown in this paper, be used to study 
these topics and more. 

Supermassive black holes, which drive AGN, are be¬ 
lieved to lie at the center of nearly every galaxy. They not 
only grow with their host galaxy, but also shape and influ¬ 
ence the galaxy and inter-galactic medium via feedback from 
their relativistic jets. So, to understand galaxies and their 


evolution we need to understand AGN evolution (eg. Fabian 
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20121. A critical key to this is identifying both young AGN 


and AGN at high-redshifts, something that is currently a 
great observational challenge. 

While AGN have been found in optical surveys out to 
z = 7.1 ( Mortlock et al.||2011 |, Ly-alpha absorption makes 
detecting them at z > 6.5 very difficult ( |Mortlock et al.| 
2011 Becker et al.|200Tl. On the other hand, the most pow¬ 


erful AGN in the radio, such as Cygnus A, are expected to 
have observed 1 GHz flux densities in the order of 15 mjy 
at z = 8, well within the capabilities of modern radio tele¬ 
scopes. Hence, not only should these sources be detectable 
in the radio, but radio also has the advantage that it is not 
obscured by dust which can hide sources at all redshifts from 
optical identification. 

While large-scale radio surveys such as the Very Large 
Array (VLA) Faint Images of the Radio Sky at Twenty- 
Centimeters (FIRST) survey ( | White et al.|[l997| ) can de¬ 
tect these sources, a method needs to be found to iden¬ 
tify them based purely on their radio properties. Based on 
an observed correlation between redsliift and spectral index 
(eg. De Breuck et al.j|200CM, searches for high-redshift radio 


galaxies tend to focus on sources with ultra-steep spectral 
indices. While this method has proven successful in finding 


sources out to z ~ 4 (eg. 

Jarvis et al. 

2001 

Cruz et al. 

2006 De Breuck et al.|2006 

), Singh et al. 

(2014 

1 found me- 


dian redshifts of only ~ 1.18 and ~ 1.57 for their two selec¬ 
tions of faint ultra-steep-spectrum (USS) sources, which is 
only slightly higher than the median redshifts of ~ 0.99 and 
~ 0.96 for the non-USS sources in their samples. Moreover, 
the physical reason for why USS sources should be at higher 


redshifts than non-USS sources remains unclear (Miley & 
Dc Brcuck 2008) and indeed recent deep observations of the 
Cosmological Evolution Survey (COSMOS) field with the 
VLA found no clear evidence that sources at higher redshift 
have steeper spectral indices |Smolcic et al. |[2014[ ). Finally, 
in a recent study using nine samples of radio sources, |Ker| 


et al. (2012) found that the fraction of 2 > 2 sources is not 


significantly higher in the sub-sample of USS sources com¬ 
pared to the full sample. Hence, finding a method with which 
high-redshift AGN can reliably be selected from radio im¬ 
ages (the focus of this paper) is very important as it would 
allow us to search for them in existing large radio surveys 
such as FIRST, current and future LOFAR (Low-Frequency 
Array; van Haarlem et al. :20l3) surveys such as the Multi¬ 
frequency Snapshot Sky Survey (MSSS) and in the future, 
the Square Kilometer Array (SKA), which might be able to 
detect these sources beyond redshift ten ( [Falcke et al.|2004 1. 

High-frequency peaked (HFP), Gigahertz peaked- 
spectrum (GPS) and compact steep-spectrum (CSS) sources 
are thought to be the youngest and smallest AGN ( |Q’Dea 
1998 Murgia et al.|[2002 Conway||2002 1. These sources are 


characterized by their peaked spectra and steep spectral in¬ 
dices (a; where a is defined as S ~ v a with S the spectral 
flux density at frequency v) above the spectral turnover. 
The HFP and GPS sources have rest-frame turnover fre¬ 
quencies above 1 GHz and linear (lobe to lobe) sizes smaller 
than ~ 1 kpc while the CSS sources typically turn over be¬ 
low 500 MHz and range in size from 1 kpc to 20 kpc. An 
evolutionary sequence is believed to exists from the HFP 


to the GPS sources and finally to the CSS sources (O’Dea 
1 1998|[Snellen et al.|2000||Tschager et al.|2003| |. Furthermore, 
it is likely that the CSS sources will evolve to eventually 


become the large FRI and FRII radio galaxies (Begelman 
1996 Snellen et al.pOOO de Vries et al.(2002 l. 


An empirical relation was found between the rest-frame 
turnover frequency u r and the linear size of nearby (z ~ 1) 


GPS and CSS source over three orders of magnitude (O’Dea 
1998 Snellen et al.pOOO). From the observed turnover fre¬ 


quency v a we can therefore estimate the linear size of a 
source by taking into account that v r = v 0 (l + z), where 
z is the redshift of the source. Using this, we can predict the 
source’s observed angular size, noting that in conventional 
cosmology (fi m = 0.3, = 0.7, Ho = 72kins -1 Mpc -1 ), a 

fixed ruler will appear to first decrease in angular size mov¬ 
ing from 2 = 0 to z ~ 1 before appearing to increase again 
beyond 2 ~ 1. We therefore expect a typical GPS source of 
100 pc, with corresponding rest-frame turnover frequency of 
2.7 GHz, to appear as a 25 mas source at 2 = 10 with an un¬ 
usually low observed turnover frequency of 245 MHz. From 
these arguments, Falcke et al. (20041 proposed an, as of yet, 
untested method to find high-redshift AGN by searching 
for compact (on scales of tens of milliarcseconds) megahertz 
peaked-spectrum (MPS; turnover frequency below 1 GHz) 
sources. 

It should be pointed out that, while our aim is to use the 
MPS method to search for high-redshift (z > 2) AGN, we 
expect that the method will also find young AGN at lower 
redshifts. For example, sources with a rest-frame turnover 
frequency of 1 GHz and corresponding size of 470 pc will have 
an observed turnover frequency of 325 MHz and angular size 
of 60 mas at 2 = 2.1. While this means that we do not expect 
all the MPS sources that are compact on a few hundred 
milliarcseconds to be at high redshifts, the sources that are 
not at high redshifts are still young AGN and are interesting 
in their own right since they help us understand how the 
AGN population evolves with redshift. 

In this paper, we publish a 325 MHz source catalogu^jof 
the Bootes field that is at least an order of magnitude better 
in terms of both resolution and noise (see Table[2j than what 
was previously available for this field with the Westerbork 
Northern Sky Survey (WENSS; Rengelink et al.||1997 l. We 
also take the first steps towards using the MPS method to 
search for high-redshift AGN by identifying 33 MPS sources 
in the image and looking at their redshift distribution. 

In Section [2] we describe the data reduction and imag¬ 
ing before describing source extraction, image quality and 
catalogue construction in Section [3] The spectral proper¬ 
ties of the sources, selecting and discussing the MPS sources 
and the affect of selecting USS sources in different frequency 
ranges on the selected sample of USS sources is presented 
in Section [4] Finally, in Section [5] we summarize our results 
and discuss future work. 


2 OBSERVATIONS AND DATA REDUCTION 

The archival data set (project code AB0976) of the Bootes 
field was taken with the VLA on 22 November 2000, 12 
December 2000, 2 January 2001 and 4 January 2001 and 
is summarized in Table [T] The observations were taken in 


1 Available in electronic format from the CDS at http://edsare.u- 
strasbg.fr/ 
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Table 1 . Summary of observations. 


Date 

Start Time 
[UTC] 

End Time 
[UTC] 

TOS a 
Bootes 
[min] 

TOS a 
3C286 
[min] 

22 Nov. 2000 

12:06:59 

20:02:56 

255 

42 

12 Dec. 2000 

12:17:30 

20:14:00 

280 

34 

2 Jan. 2001 

10:22:56 

18:51:20 

285 

38 

4 Jan. 2001 

10:47:40 

18:43:33 

381 

41 


Notes: a The time on source (TOS). 


the P-band with the VLA in A configuration, using 27 an¬ 
tennae on every day except 4 January 2001 when only 26 
antennae were available. The observations covered a single 
pointing in the Bootes held centered at RA 14:32:05.72 and 
DEC +34:16:47.5 using 3.3 s integrations with a total in¬ 
tegration time on the held of 20 hours spread over the four 
days. The data were taken in spectral line mode with two in¬ 
termediate frequencies (IFs) at 327.53 MHz and 321.53 MHz 
each with two polarizations (LL and RR). We will take the 
average frequency of 324.525 MHz of the two IFs as the rep¬ 
resentative frequency. Each IF was split into 16 channels of 
394 KHz each, giving a bandwidth of 6.3 MHz per IF and 
total bandwidth of 12.6 MHz. Observations of the calibrator 
source 3C286, used for phase, amplitude and bandpass cal¬ 
ibration, were irregularly interleaved between those of the 
target held. Typically, 3C286 was observed for 2.5 minutes 
every 25 minutes. 

The initial data reduction was done using the AIPS 
(Greisen 19901 software package. Data were hagged using 
both the automatic bagging task rflag and the manual 
flagging task tvflg before and after calibration. After cal¬ 
ibration, the calibrated data were imported into CASi0 for 
self-calibration, imaging and primary beam correction. Five 
rounds of phase only self-calibration was done, in which 
the solution interval was decreased from 20 to 1 min. Dur¬ 
ing each round, a clean component model of the brightest 
sources in the held was derived from the previous best im¬ 
age and used as an input model. The final image of the 
held (the inner 1.2° of which is shown in Fig. [l]) was made 
using 1.5arcsec pixels, Brigg’s weighting with a robustness 
parameter of 0.0 and 128 w-projection planes. While the 
VLA P-band primary beam has a full-width half-maximum 
(FWHM) of 2.5°, the hnal image produced using CASA’S 
CLEAN task was made with a diameter of 4.1° to increase the 
number of sources in the image. The image has a resolution 
of 5.6 x 5.1 arcsec and central noise of ~ 0.2 mjy beam -1 in¬ 
creasing to ~ 0.8 mjy beam -1 at the edge of the image after 
the hux correction has been applied (see Section 3.31. 


3 IMAGE QUALITY AND CATALOGUE 
CONSTRUCTION 

3.1 Source Extraction 

Source extraction was done using the pybdsm source de¬ 
tection package]^] pybdsm works by identifying all pixels 


2 http://casa.nrao.edu 

3 http://dl.dropboxusercontent.com/u/1948170/html/index.html 


in the image that are above the pixel detection threshold 
and adding all of the surrounding pixels that are above the 
island detection threshold together to form islands of emis¬ 
sion. Multiple Gaussians are then fitted to each island before 
grouping the Gaussians together into individual sources. We 
used the default values of three and five times the the local 
root mean square (mis) noise for the island and pixel thresh¬ 
olds respectively, where the local rms noise was calculated 
in a sliding box of 50 x 50 pixels that was moved across the 
image in steps of ten pixels. 

For each source, the flux density reported is the sum 
of the flux densities of the Gaussians in the source, the flux 
density errors are the 1 <t statistical errors calculated using 
the formula in Condon (1997) and the source position is the 
the centroid of the Gaussians. In total 1377 sources made 
up of 1586 Gaussians were recovered from the image. Af¬ 
ter visual inspection, seven sources were rejected because 
they were found to be artefacts near bright sources. Of the 
remaining 1370 sources, 1251 were fit by a single Gaussian. 


3.2 Source Matching 


To determine the spectral properties of the sources and 
check the flux density scale of our image (hereafter referred 
to as the VLA-P image), all of the sources were matched to 
the 1.4 GHz FIRST, 1.4 GHz National Radio Astronomy Ob¬ 


servatory (NRAO) VLA Sky Survey (NVSS; Condon et al. 
|1998[ ), 1.38GHz ( |de Vries et al. |[2002[ hereafter referred to 
as the deVries catalogue), 325MHz WENSS and 153MHz 
I Will iams et al.| [2013| hereafter referred to as the Williams 
catalogue) catalogues. A summary of the catalogues (includ¬ 
ing the VLA-P catalogue for comparison) is given in Table[2] 
Matching was done with the stilts ( Taylor| [2006) software 
package using the task tskymatch2. Using the central po¬ 
sition of each sources in the VLA-P catalogue and a search 
radius of half the FWHM of the catalogue with the lowest 
resolution, tskymatch2 returned all sources in the matched 
catalogue separated by less than the search radius from the 
input sources. 

Since our aim is to search for high-redshift AGN, which 
are expected to have linear sizes of at most a few hundred 
milliarcseconds, we matched the VLA-P sources to the cat¬ 
alogue with the highest resolution at each frequency. How¬ 
ever, while FIRST has the highest resolution at 1.4 GHz, 562 
sources could not be matched to a FIRST source because 
they were too faint to be detected in FIRST due to their 
steep spectral indices. Consequently, the VLA-P sources 
were also matched to the deVries catalogue, which, while 
having a lower resolution than FIRST (see Table [2|, also 
has a significantly lower rms noise. Hence, 417 sources that 
could not be matched to a FIRST source could be matched 
to a source in the deVries catalogue. 

In order to exclude confused sources and sources that 
are resolved in one of the catalogues, only one-to-one 
matches were accepted. In other words, if a source in the 
VLA-P catalogue was matched to two or more sources in 
another catalogue, or vise versa, the matches were rejected. 
To exclude resolved sources, VLA-P sources with a major 
or minor axis greater than 10 arcsec were excluded from all 
further analysis. While the beam size in the VLA-P image 
is ~ 6 arcsec, a value of 10 arcsec was used to accommo¬ 
date the errors on the fitted major and minor axes of the 
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BOOTES 



Figure 1. Inner 1.2° of the VLA-P image. The full image has a radius of 2.05°, a resolution of 5.6 X 5.1arcsec and central noise of 
~ 0.2 mjy beam -1 . 


Table 2. Catalogues matched to the VLA-P catalogue. 


Catalogue 

Name 

Frequency 

[MHz] 

Resolution 

["] 

rms Noise a 
[mJy beam -1 ] 

FIRST 

1400 

5.4 x 5.4 

0.15 

NVSS 

1400 

45 x 45 

0.45 

deVries 

1380 

13 x 27 

0.03 

VLA-P 

325 

5.6 x 5.1 

0.2* 

WENSS 

325 

54 x 54 

3.6 

Williams 

153 

25 x 25 

2.0* 


Notes: a The values quoted are the typical catalogue noise except 
for the two marked with a '*’, which are the noise at the center 
of the image. 


faint sources. Similarly, VLA-P sources that were matched 
to a source with a major or minor axis greater than 7 arcsec 
in FIRST or matched to a sources that is flagged as being 
extended in WENSS were excluded. In total 283 resolved 
sources were excluded. Matching to multiple catalogues at 


the same frequency (FIRST, deVries and NVSS), we exclude 
323 variable sources by removing sources that showed a flux 
density difference of more than 20 per cent in any two of 
the catalogues. Since NVSS’s resolution is about eight times 
lower than that of FIRST, this also ensures that none of 
the sources have extended structure that is resolved out in 
one of the higher resolution catalogues. Finally, all of the 
matches were visually inspected to check for potential in¬ 
correct matches. This resulted in 18 matches to 13 VLA-P 
sources being rejected. 

3.3 Flux Density Scale 

Since the central frequencies of both the VLA P-band and 
WENSS catalogues are 325 MHz, we matched the two cata¬ 
logues so that we could check the absolute flux density scale 
and the primary beam correction. When doing this, only 
single isolated non-extended sources were used. Addition¬ 
ally, sources with a signal-to-noise ratio (SNR; defined as 
the peak flux density divided by the island rms) below ten 
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Distance from the phase center (degrees) 


Figure 2. The uncorrected flux density ratio (Sy LA—P *^wenss) 
plotted as a function of the distance from the phase center for the 
sources used to check the flux density scale. The dashed line shows 
the function y = 1 while the solid line is the function fitted to the 
data to correct for the VLA-P flux density offset. The function 
has a value of 1.0 at the phase center, increasing to 3.72 at the 
edge of the image. 

in either the VLA-P or WENSS catalogue were excluded 
from the analysis. 

To check the flux density scale, the flux density ratio 
of the VLA-P with WENSS c! VLA ~ P ) was plotted for each 
pair of matched sources as a function of the source SNR and 
the angle of the source with respect to the phase center in 
order to check for a potential irregular beam shape. In addi¬ 
tion, the VLA-P flux density was also plotted as a function 
of the WENSS flux density for each pair of matched source. 
No trends were visible in any of these plots. Finally, the flux 
density ratio was plotted as a function of the distance of 
the source from the VLA-P phase center (shown in Fig. [2|. 
From the figure it is clear that sources further from the phase 
center (specifically those beyond the primary beam) had a 
larger flux density ratio, indicating that the flux densities 
at the edge of the VLA-P image were over estimated. This 
problem was likely caused by an error in the model of the 
VLA P-band primary beam that CASA used to do the pri¬ 
mary beam correction. 

To correct for this offset, an exponential function of the 
form 

- = k 1 +k 2 e 3 (1) 

O WENSS 

was fitted to the data using a non-linear least squares fitting 
routine (shown as the solid line in Fig.[2|, where each point 
was weighted by its error. In the function, ki, k 2 and k 2 are 
constants with best-fit values of 1.00±0.03, 0.005±0.002 and 
2.93 ± 0.19 respectively and r is the distance of the source 
from the phase center in degrees. The fitted function has a 
reduced chi-squared of 1.34, a value of 1.003 at r = 0° and 
1.22 at r = 1.25°, the FWHM of the VLA primary beam, 
increasing to 3.72 at the edge of the image. This shows that 
the flux density offset was only a significant problem outside 
the primary beam. Using Equation [l] and the condition 

Ocorrected_VLA —P 7 ..., 

-A- = « 1, (2) 



Figure 3. Log-log plot of the corrected VLA-P flux density as a 
function of the WENSS flux density for the sources used to check 
the flux density scale. The line y = x was added as a guide. 

all of the flux densities, flux density errors and rms values, 
were corrected using their distance from the phase center 
and the errors of the fitted parameters. 

Fig-0 shows a plot of the corrected VLA-P flux den¬ 
sities as a function of the WENSS flux densities for the 
sources used to check the flux density scale. As can be seen, 
the VLA-P flux densities agree well with those of WENSS. 
Since the mean of the corrected flux density ratios is 1.05, 
we conclude that no systematic offset between the VLA-P 
and WENSS flux densities remained after the flux density 
correction. 

3.4 Catalogue Construction 

A sample of the final catalogue is given in Table 0 the 
full catalogue containing 1370 sources is available from the 
CDf-0 In the catalogue, the flux densities, peak flux den¬ 
sities and rms values have been corrected as described in 
Section [3.3| The columns in the table are: (1) Source name, 
(2,3) right ascension (RA) and uncertainty, (4,5) declina¬ 
tion (DEC) and uncertainty, (6) integrated 342.5 MHz flux 
density and uncertainty, (7) peak 342.5 MHz intensity and 
uncertainty, (8,9) deconvolved major- and minor-axis, (10) 
position angle, (11) local rms noise, (12) number of Gaus- 
sians out of which the source is composed. In all cases, values 
that could not be calculated are indicated as For sources 
composed of more than one Gaussian, the individual Gaus- 
sians are named alphabetically and given below the source. 

4 RESULTS AND DISCUSSION 
4.1 Spectral Properties 

Two spectral indices, ai ow and ahigh, were calculated for 
each source. The low frequency spectral index, ai ow , was 
calculated between the Williams (153 MHz) and VLA-P 
(325 MHz) catalogues. The high frequency spectral index, 

4 http://cdsarc.u-strasbg.fr/ 
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Table 3. Example entries of the online catalogue. 


Source ID 

RA 

CRA 

DEC 

ODEC 

Si 

S p 

a 

b 

<t> 

rms 

-^Gauss 


[deg] 

["] 

[deg] 

["] 

[mjy] 

[mjy beam -1 ] 

["] 

["] 

[deg] 

[mjy beam -1 ] 


(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(id 

(12) 

J143821+340235 

219.5863 

0.6 

34.0432 

0.6 

3.6±0.9 

2.7±0.6 

8.2±2.0 

4.7±0.7 

137±18 

0.5 

1 

J143813+335310 

219.5529 

0.7 

33.8861 

0.4 

6.6±1.0 

3.4±0.6 

9.1±1.6 

6.1±0.8 

108±21 

0.5 

1 

J143622+335939 

219.0909 

0.1 

33.9942 

0.1 

223.5±11.0 

70.5±3.5 

14.2±0.1 

7.1±0.1 

9±1 

0.4 

3 

J143622+335939a 

219.0911 

0.1 

33.9955 

0.1 

107.0±5.3 

32.7±1.7 

11.5±0.1 

8.1±0.1 

16±2 

0.4 

- 

J143622+335939b 

219.0906 

0.1 

33.9926 

0.1 

79.9±4.0 

63.0±3.1 

6.7±0.1 

5.4±0.1 

117±2 

0.4 

- 

J143622+335939c 

219.0908 

0.1 

33.9940 

0.1 

36.6±1.9 

22.4±1.2 

7.8±0.1 

6.0±0.1 

127±4 

0.4 

- 

J142208+342122 

215.5341 

0.7 

34.3561 

0.4 

10.7±4.3 

6.3±2.6 

7.9±1.7 

6.0±1.0 

87±31 

1.1 

1 

J143948+330413 

219.9511 

0.8 

33.0702 

0.4 

26.1±9.2 

5.8±2.2 

8.1±1.9 

5.5±1.0 

172±26 

1.0 

2 

J143948+330413a 

219.9516 

0.7 

33.0690 

0.5 

8.0±3.2 

5.0±2.0 

7.7±1.8 

5.9±1.1 

69±38 

1.0 

- 

J143948+330413b 

219.9509 

0.8 

33.0709 

1.0 

18.1±6.4 

5.3±2.1 

12.8±2.7 

7.6±1.3 

38±22 

1.0 

- 


cthigh, was calculated between the VLA-P (325 MHz) and 
FIRST (1400 MHz) catalogues unless the source could not 
be matched to a source in FIRST (for the reason explained in 
Section [±2] but could be matched to a source in the deVries 
catalogue. In that case cthigh was calculated between the flux 
densities of the VLA-P (325 MHz) and deVries (1380 MHz) 
catalogues. 

The anigh values calculated using FIRST were found to 
have a median value of -0.72 and interquartile range from 
—0.91 to —0.52. These values are consistent with those found 
by |Smolcic et al.| ( |2014[ ) and |Owen et al.| ( |2009| ) who made 
VLA P-band images of the VLA COSMOS and Spitzer 
Wide-area InfraRed Extragalactic Survey (SWIRE) fields 
with similar sensitivity and resolutions as our image. This 
serves to further confirm our flux density scale. The me¬ 
dian value and interquartile range for aiow was found to be 
-0.56 and from -0.83 to -0.30. Considering that cthigh and 
aiow have mean values of —0.685 ±0.007 and —0.505 ±0.033 
respectively, is is clear that the spectral index flattens to¬ 
ward lower frequencies, as was reported by among others 


Williams et al. (2013) and van Weeren et al. (20141. This 


is possibly the result of spectral aging which results in the 
sources having a steeper spectrum a higher frequencies. 


4.2 MPS Sources 

To select the MPS sources, a colour-colour diagram (Fig. [4j) 
was made by plotting cthigh against ai ow for the 198 sources 
which could be matched to both the Williams and either the 
FIRST or the deVries catalogues. In the figure, the sources 
in the bottom left and top right quadrants have, respec¬ 
tively, negative and positive spectral indices across the en¬ 
tire frequency range from 153 MHz to 1.4 GHz and contain 
85 and 3 per cent of the sources respectively. The remaining 
12 per cent of the sources lie in the bottom right quadrant 
and have convex (peaked) spectra. In order avoid confusion, 
error bars are not shown in the figure. As a guide it is noted 
that the median errors on the points are 0.1 and 0.43 for 
cthigh and aiow respectively. It should be noted that from 
the size of the median errors, it is clear that the higher er¬ 
rors on ai ow are caused by higher flux density errors in the 
Williams catalogue compared to those in FIRST. 

Since the MPS sources have a spectral peak below 
1 GHz, we selected these sources to have ahigh < —0.5 and 
ahigh < 1-5 atom — 0.5 (the region below and to the right of 


the dotted line in Fig. [4|. The selection criteria for ai ow 
was used in order to not only select the sources with a 
clear spectral peak, but also select the sources that flat¬ 
ten towards lower frequencies, which could be indicative of 
a turnover frequency below 153 MHz. While the selection 
criteria ahigh < l-5aiow — 0.5 is somewhat arbitrary, it does 
ensure that ahigh is at least one and a half times ai ou , for 
the sources that show a flattening in their spectra. In addi¬ 
tion, this also removes the effect of the spectral index flat¬ 
tening towards lower frequencies (see Section 4.1l from the 
selection. In total, 33 MPS sources were selected and pre¬ 
sented in Table [4] We note that when fitting Gaussians to 
the sources, pybdsm fit a single Gaussian to each of the 33 
MPS sources. In the table, the columns correspond in name 
to the columns in Table [3] except for the following columns: 
(8) low frequency spectral index, (9) high frequency spec¬ 
tral index, (10) redshift calculated using the EAZY code, 
(11) redshift calculated using the LTR code (a description 
of how columns 10 and 11 were calculated are given be¬ 
low). While the low frequency spectral index was calculated 
between 153 MHz and 325 MHz, the high frequency spectral 
index was calculated between 325 MHz and 1400 MHz unless 
the value is marked with a in which case it was calcu¬ 
lated between 325 MHz and 1380 MHz. Of the 33 sources in 
Table [4] 18 have spectra that flatten ( ai ow ^ 0) while the 
remaining 15 sources have peaked spectra. Due to the errors 
on the spectral indices, only two of the 15 sources might 
unambigously be consistent with a peaked spectrum. 

Photometric redsliifts ( z p hot ) were determined for the 
MPS sources by identifying candidate counterparts in the 


NOAO Deep Wide-Field Survey (NDWFS; Jannuzi & Dey 
1999]) I-band images using either the FIRST position or, if 
the source could not be matched to a FIRST source, the 
VLA-P position. We then matched the I-band sources to, 
and extracted flux density measurements from, the point 


spread function matched photometric catalogues (Brown 


et al. 2007) which comprise fluxes from the following sur¬ 


veys: NDWFS (B_W, R, I and K bands), the Flamingos Ex¬ 


tragalactic Survey (FLAMEX, J and Ks bands; Elston et al. 
2006), the zBootes survey (z’ band; |Coo l [200 7), the Spitzer 


Deep Wide Field Survey (SDWFS, 3.6, 4.5, 5.8 and 8.0 mi- 
Ashby et al. 2009), the Galaxy Evolution Ex¬ 


cron images; 


plorer GR5 survey (GALEX, near- and far-ultraviolet bands; 
|Morrissey et al.|2007[) and the MIPS AGN and Galaxy Evo¬ 


lution Survey (MAGES, 24 micron image; Jannuzi et al. 
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Figure 4. Colour-colour diagram for sources with a one-to-one match to the Williams and either the FIRST or the deVries catalogues. 
In the plot, aq ow was calculated between 153 MHz and 325 MHz while, was calculated between 342.5 MHz and 1400 MHz for the 

sources marked as filled circles and between 342.5 MHz and 1380 MHz for sources marked as crosses. The dashed lines indicate spectral 
indices of zero and the MPS sources were selected from the region below and to the right of the dotted line. While error bars are not 
shown to avoid confusion, the median error bar size is illustrated at the bottom right of the figure. 


2010). Finally, the resulting spectral energy distribution 


(SED) was fitted for z p hot using the LRT code from As- 
sef et al. (20081 and EAZY code ( Brammer et al.||2008 ) for 


comparison. A more detailed description of the matching 
and SED fitting of the sources in the field will be given in 
Williams et al. (in preparation). 


In total, we were able to determine redshifts (given in 
Table [4| for 24 of the 33 sources. Of the nine sources with¬ 
out redshifts, five lie outside the multiwavelength coverage 
while the remaining four were too faint to be matched. 
Except for the two confused sources J143359+340421 and 
J143542+341321, each of which have two optical matches, 
and J142917+332626, the values output by the LRT and 
EAZY codes agree well with each other and range between 
0.14 and 3.16. While the LRT code, which includes an empir¬ 
ical AGN SED template in the fitting and therefore fits AGN 
spectra better, gives a mean and median redshift of 1.3 and 
1.0 respectively, the mean and median redshift output by the 
EAZY code is 1.4 and 1.1. On average, both the LRT and 
EAZY codes agree that the redshift of the flattening sources 


are 0.1 units higher than the peaked sources. Similarly, the 
median redshift of the flattening sources are 0.4 units higher 
than that of the peaked sources for the EAZY code, and 0.3 
units higher for the LRT code. While the median redshift 
error of the EAZY code is 0.1 units for both the flattening 
and peaked sources, and could therefore explain the differ¬ 
ence between the average redshifts, this does add weight to 
the argument that the flattening sources could peak towards 
lower frequencies and could therefore be at higher redshifts 
than the peaked sources. 


Since our current selection allows sources with an ob¬ 
served turnover frequency of ~ 325 MHz to be selected, and 
these sources are expected to have a linear and angular size 
of ~ 2.3 kpc and ~ 1.3arcsec respectively at z = 0.1, it 
is not surprising that our selection contains sources at low 
(z < 1) redshifts. It should however be pointed out that 
these sources are classic CSS sources and hence likely to be 
young AGN ( Q’Dea|1998 Murgia et al.|2002 |. Our expecta¬ 
tion is that the high-redshift sources should all be compact 
on scales of tens of milliarcseconds, which is something that 























8 Rocco Coppejans et al. 


Table 4. List of MPS sources. 


Source ID 

RA 

CRA a 

DEC 

^DEC a 

Si 

rms 

oq 

ow 

b 

CH-high 

^EAZY d 

~ d 

ZLRT 


[deg] 

["] 

[deg] 

H 

[mjy] 

[ mJy beam -1 ] 







(i) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 


(9) 

( 10 ) 

( 11 ) 

J143902+351652 

219.7583 

0.6 

35.2812 

0.5 

8.2 ± 2.2 

0.6 

0.21 

± 

0.84 

-1.16 ± 0.19 

- 

- 

J143833+331845 

219.6358 

0.2 

33.3126 

0.2 

13.9 ± 3.0 

0.5 

0.44 

± 

0.61 

-0.54 ± 0.15 

- 

- 

J143813+335310 

219.5529 

0.7 

33.8861 

0.4 

6.5 ± 1.0 

0.5 

-0.42 

± 

0.63 

-1.19 ± 0.12* 

i 'tqo+0.264 
t. i • 7 °_o_269 

2.06 

J143657+360144 

219.2376 

0.5 

36.0289 

0.6 

9.6 ± 3.8 

0.9 

-0.40 

± 

0.67 

-1.12 ± 0.27 

- 

- 

J143628+345219 

219.1147 

0.4 

34.8719 

0.5 

4.7 ± 0.5 

0.3 

-0.07 

± 

0.55 

-0.74 ± 0.08* 

0 Q5Q+ U Uy& 

U.OOO_q 0 g 4 

0.77 

J143624+355419 

219.0983 

0.2 

35.9052 

0.5 

12.7 ± 3.9 

0.7 

0.07 

± 

0.60 

-0.55 ± 0.21 

- 

- 

J143619+350156 

219.0805 

0.1 

35.0321 

0.1 

16.2 ± 1.3 

0.3 

0.29 

± 

0.42 

-0.71 ± 0.07 

0 899 +0 ' 051 
u - oyy _ 0 . 051 

0.74 

J143542+330225 

218.9244 

0.4 

33.0402 

0.4 

5.1 ± 1.0 

0.4 

-0.41 

± 

0.62 

-1.15 ± 0.14* 

9 

z.ouo_o 165 

2.10 

J143542+341321 

218.9259 

0.1 

34.2224 

0.1 

5.6 ± 0.4 

0.2 

0.44 

± 

0.67 

-0.85 ± 0.06 

2.850+°'J°® or 

5.30 or 












2-856™ 

5.09 

J143535+325959 

218.8958 

0.3 

32.9998 

0.4 

6.9 ± 1.2 

0.4 

-0.03 

± 

0.62 

-0.99 ± 0.12 

3 162 -1-0 ' 191 
o.ioZq 156 

3.04 

J143505+334721 

218.7707 

0.0 

33.7893 

0.0 

17.1 ± 0.9 

0.2 

-0.14 

± 

0.29 

-0.72 ± 0.05 


- 

J143359+340421 

218.4956 

0.1 

34.0725 

0.1 

4.2 ± 0.4 

0.2 

-0.23 

± 

0.54 

-0.89 ± 0.07* 

3.615+°'J°i or 

0.62 or 












A - ool _0.218 

1.07 

J143339+354013 

218.4135 

0.1 

35.6702 

0.2 

10.5 ± 1.5 

0.4 

0.07 

± 

0.44 

-0.84 ± 0.11 

0 631 +U,U69 

U.Ool _ 0 07 1 

0.54 

J143330+355042 

218.3741 

0.1 

35.8449 

0.2 

14.6 ± 2.9 

0.5 

0.16 

± 

0.48 

-0.51 ± 0.14 

o 09 - 1 -|-2.032 
z.ozi — i 536 

1.37 

J143230+353641 

218.1256 

0.1 

35.6113 

0.3 

9.0 ± 1.1 

0.3 

-0.14 

± 

0.43 

-0.87 ± 0.09 

9 01 rj+0.253 

^.otu_o 248 

2.36 

J143223+324940 

218.0949 

0.0 

32.8277 

0.1 

27.3 ± 3.9 

0.4 

0.02 

± 

0.35 

-0.62 ± 0.10 

- 

- 

J143215+350608 

218.0621 

0.2 

35.1023 

0.2 

4.0 ± 0.4 

0.2 

-0.30 

± 

0.55 

-1.26 ± 0.08* 

0 747 + 0 - 043 
u '' ^' —0.046 

0.86 

J143138+353215 

217.9064 

0.3 

35.5375 

0.3 

5.5 ± 0.7 

0.3 

- 0.11 

± 

0.69 

-0.96 ± 0.10 

H OO7+0-091 
_0.069 

0.69 

J143123+331626 

217.8472 

0.0 

33.2738 

0.0 

51.5 ± 3.1 

0.2 

-0.13 

± 

0.29 

-0.81 ± 0.05 

- 

- 

J143051+342614 

217.7120 

0.0 

34.4372 

0.0 

29.0 ± 1.3 

0.2 

0.49 

± 

0.34 

-0.61 ± 0.05 

2 364 +ut)iib 

^■ OU4 _0.536 

2.98 

J143001+340746 

217.5049 

0.1 

34.1295 

0.1 

7.4 ± 0.5 

0.2 

0.52 

± 

0.80 

-0.78 ± 0.05 

U.lOO _ 0 05 5 

0.19 

J142938+343903 

217.4085 

0.0 

34.6507 

0.0 

14.5 ± 0.8 

0.2 

-0.36 

± 

0.36 

-1.10 ± 0.05 

0 139+ 0037 

U. -LO«y— 0.037 

0.15 

J142941+330552 

217.4195 

0.1 

33.0978 

0.2 

8.6 ± 1.0 

0.3 

0.27 

± 

0.59 

-0.76 ± 0.09 

1 71 q+0-335 
t. i xy _o.337 

1.38 

J142917+332626 

217.3226 

0.0 

33.4407 

0.0 

17.3 ± 1.1 

0.3 

0.60 

± 

0.37 

-0.64 ± 0.06 

1 kqo+0.322 
t. 000 _q.290 

2.49 

J142905+354425 

217.2691 

0.1 

35.7402 

0.2 

13.0 ± 2.5 

0.5 

- 0.00 

± 

0.48 

-0.59 ± 0.14 

0 809 +O '° 84 

u.ouy_ 0 .081 

0.84 

J142906+343326 

217.2764 

0.2 

34.5571 

0.1 

7.9 ± 0.6 

0.3 

-0.17 

± 

0.45 

-0.77 ± 0.06 

n * 377 + 0.019 
u -°‘ 1 -0.032 

0.33 

J142836+353154 

217.1479 

0.1 

35.5315 

0.1 

29.7 ± 4.2 

0.4 

- 0.10 

± 

0.35 

-0.81 ± 0.10 

2 190”^°' 254 
Zr.xyu_o 235 

2.21 

J142719+352324 

216.8305 

0.3 

35.3900 

0.2 

5.1 ± 1.0 

0.4 

-0.61 

± 

0.57 

-1.72 ± 0.17* 

1.69012: (224 

1.56 

J142717+351756 

216.8197 

0.2 

35.2989 

0.2 

7.7 ± 1.2 

0.4 

0.35 

± 

0.70 

-0.67 ± 0.11 

0.820l":““ 

0.87 

J142722+333123 

216.8406 

0.2 

33.5230 

0.2 

5.2 ± 0.7 

0.3 

0.04 

± 

0.82 

-0.76 ± 0.10 

- 

- 

J142608+341658 

216.5339 

0.3 

34.2829 

0.2 

5.5 ± 0.7 

0.3 

-0.11 

± 

0.69 

-0.67 ± 0.09 

1 329 +0 104 

i.ozy_Q 104 

1.13 

J142602+354639 

216.5065 

0.5 

35.7775 

0.4 

11.1 ± 4.0 

0.8 

0.13 

± 

0.72 

-0.51 ± 0.25 

U.oo9 _ 0 Q40 

0.68 

J142310+333033 

215.7924 

0.4 

33.5092 

0.2 

18.5 ± 6.8 

0.9 

-0.04 

± 

0.60 

-0.71 ± 0.26 

- 

- 


Notes: a RA and DEC errors smaller than 0.04arcsec are rounded to 0.0. b oq ow was calculated between 153 MHz and 325 MHz. c If 
cthigh was calculated between 325 MHz and 1380 MHz, the value is marked with a otherwise cthigh was calculated between 342.5 MHz 
and 1400 MHz. d The entries containing two values are because the radio source was matched to two optical sources. In the table, the 
redshift of both optical matches are given. 


we will test in our upcoming paper Coppejans et al. (in 
preparation), and appear as quasars at optical wavelengths 
( |Q’Dea| 1998| |1990) . What is encouraging is that both codes 
agree that five of the sources for which we have redshifts are 
at 2 > 2. Furthermore, given the correlation between opti¬ 
cal magnitude and redshift for the host galaxies of sources 
with radio jets (O’Dea 1998 Rocca-Volmerange et al.|2004 


Miley fe De Breuck||2008 1 , we expect that the four sources 

for which we could not find redshifts because they are too 
faint to be matched are likely at higher redshifts than the 
sources with known redshifts. 


4.3 USS Sources 


USS source have been selected using different spectr al cuts 


and frequency ranges eg., < 


-1.1 (Wieringa & 


ICatgert 1992), al f, 5 ,™ < —0.981 (Blundell et al. 

1998), 

0 * 843 MH 7 —1.3 (De Breuck et al. 2004), ot\ 5 ^ “C —l.C 

(Cruz et al. 2006), aJosMHz < -1.0 

Broderick et al. 

2007) 

and £* 325 mhz < -10 (Singh et al. 

2014). Defining USS 


sources as those with a spectral index steeper than -1.3, 
we find that only 4.2 per cent of the sources for which cthigh. 
could be calculated can be classified as USS sources. Simi¬ 
larly 4.9 per cent of the sources for which ai ow could be cal¬ 
culated can be classified as USS sources. This is once again 


consistent with both Smolcic et al. (2014) and Owen et al. 


12009) who do not find a large population of USS sources in 
their samples. 


Looking only at the sources in the colour-colour dia¬ 
gram (for which both cthigh and ai ow are available), we find 
that eight sources can be classified as USS based on their 
value for OLhigh and another seven based on their value of 
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otiow What is interesting is that only three sources have 
both ahigh and aio w smaller than -1.3 and will appear in 
both selections. Hence, if we make two selections of USS 
sources, one based on ahigh and the other based on ai ow , 
and take into account the error bars on the spectral indices, 
only 25^8 per cent of the sources in the combined selection 
appear in both of the individual selections. Changing the se¬ 
lection criteria for USS sources from a spectrum steeper than 
-1.3 to a spectrum steeper than -1.0, changes the number to 
19^g 2 percent. To check if these low values are because ai ow 
is on average 0.2 units higher than ahigh (see Section 4.11, 
we made new selections use a cut that is 0.2 units higher for 
aiow than that of ahigh- For both cases, the percentage of 
sources that appear in both selections were between two and 
five per cent lower than the values given above, indicating 
that the low percentages are not caused by spectral flatten¬ 
ing at low frequencies. Hence, in our case, even when using 
the same spectral cut, selecting USS sources in different fre¬ 
quency ranges does not select the same group of sources and 
care should be taken when comparing different selections of 
USS sources. 


for searching for MPS sources. As pointed out in Section 
|4.2| the large errors on the low frequency spectral indices 
are the result of large flux density errors in the 153 MHz 
Williams catalogue that was made using observations by 
the Giant Metrewave Radio Telescope (GMRT). Since it is 
expected that LOFAR will produce images with lower noise 
and higher resolution than what can be achieved with the 
GMRT at these frequencies, it will significantly improve the 
reliability with which MPS sources can be selected. 

Finally, in Section |4.3| we found that selecting USS 
sources in different frequency ranges using the same spec¬ 
tral cut does not select the same group of sources and care 
should be taken when comparing different selections of USS 
sources. 
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5 SUMMARY AND FUTURE PROSPECTS REFERENCES 


In this paper, we presented a 325 MHz VLA image of the 
NOAO Bootes field. Both the resolution and the noise of 
image are an order of magnitude better than what was pre¬ 
viously available for the field at this frequency. Matching our 
sources to the FIRST, |de Vries et al.| ( |2002| l and |Williams| 


et al. (2013) catalogues of the field we calculated spectral 


indices for the sources between 153 MHz and 325 MHz and 
between 325 MHz and 1.4 GHz. We then used the spectral 
indices to make a radio colour-colour diagram of the uncon¬ 
fused point sources. 

In Section [l] we described an untested method of find¬ 
ing young and high-redshift AGN by selecting sources with a 
peaked spectra in the megahertz frequency range, the mega¬ 
hertz peaked-spectrum (MPS) sources. Using our colour- 
colour diagram, we identified 33 MPS sources of which 15 
have peaked and 18 have flattening spectra. Of the 33 MPS 
sources, we were able to determine redshift values for 24 
ranging between 0.1 and 3.2. Considering that five of the 
sources are at 2 > 2, that we expect the sources that are 
compact on scales of tens of milliarcseconds to be at the 
highest redshifts while we could only select sources that are 
compact on a scale of ~ 5 arcsec and that the four sources 
which were too faint to be matched are likely at a > 3, 
there is encouraging evidence that the MPS method can be 
used to search for high-redshift sources. In the second pa¬ 
per in this series, Coppejans et al. (in preparation), we will 
take the next step in testing the MPS method by presenting 
observations of the sources with the European Very-long- 
baseline interferometry Network (EVN) to determine their 
angular sizes and combine that with our redshift estimates 
to test whether the sources at higher redshifts are compact 
on scales of tens of milliarcseconds as we expect. 

At the moment, the greatest challenge to finding MPS 
sources is the lack of high resolution (< 10 arcsec) low noise 
(<1 mjy beam -1 ) radio maps below 1 GHz. Since LOFAR 
will be able to make such maps in the near future and covers 
the frequency range between 30 and 240 MHz, it will be ideal 
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